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ABSTRACT 

The cosmological dark matter field is not completely described by its hierarchy of Ap¬ 
point functions, a non-perturbative effect with the consequence that only part of the 
theory can be probed with the hierarchy. We give here an exact characterization of the 
joint information of the hierarchy within the lognormal field. The lognormal field is the 
archetypal example of a field where this effect occurs, and, at the same time, one of the 
few tractable and insightful available models to specify fully the statistical properties 
of the evolved matter density field beyond the perturbative regime. Nonlinear growth 
in the Universe in that model is set letting the log-density field probability density 
functional evolve keeping its Gaussian shape, according to the diffusion equation in 
Euclidean space. We show that the hierarchy probes a different evolution equation, 
the diffusion equation defined not in Euclidean space but on the compact torus, with 
uniformity as the long-term solution. The extraction of the hierarchy of correlators 
can be recast in the form of a nonlinear transformation applied to the held, ’wrap¬ 
ping’, undergoing a sharp transition towards complete disorder in the deeply nonlinear 
regime, where all memory of the initial conditions is lost. 
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1 OVERVIEW 

According to the successful cosmological A cold dark mat¬ 
ter (ACDM) model, the evolution of the matter density field 
can be represented by the Vlasov-Poisson system in an ex¬ 
panding Universe, starting from Gaussian initial conditions. 
The study of this nonlinear system is a complex endeavour, 
typically addressed with perturbation theory (PT) in vari¬ 
ous forms or V-body numerical simulations nowadays, mo¬ 
tivations being not only better understanding of structure 
formation, but also the fact that most cosmological probes 
of the low redshift Universe and of its accelerated expansion 
(Albrecht et al. 2006; Weinberg et al. 2013) must rely, in 
order to fulfil their promises, on a thorough qualitative and 
quantitative understanding of the field’s statistical proper¬ 
ties. 

Complete specification of a field statistics requires a 
probability density ’functional’, a self-consistent prescrip¬ 
tion assigning probabilities to joint occurrences at all points 
(or, for all practical purposes, at the different resolution ele¬ 
ments, or ’pixels’). In the cosmological framework this func¬ 
tional is given formally by an intractable Liouville equation, 
describing conservation of probability along trajectories in 
phase space. Perturbative methods such as PT can trace suc¬ 
cessfully some simple properties of that functional close to 


the linear regime, such as two-point statistics, while A-body 
numerical methods provide single field realizations only. In¬ 
sights gained from greatly simplified models can thus often 
prove invaluable. The lognormal density field (Coles & Jones 
1991) is one of these models. It is defined as the exponential 
of a Gaussian random field. The Gaussian log-density field is 
fully specified by requesting the resulting density two-point 
function to match that of the cosmological model. 

One wonders how a parameter-free, local transforma¬ 
tion of a Gaussian field can possibly explain many qual¬ 
itative properties of the evolving matter field. The most 
compelling argument at the time of writing is probably 
that this field exhibits hierarchical form (Fry 1984; Pee¬ 
bles 1980) in the weakly nonlinear regime of structure for¬ 
mation, and reduces naturally to a Gaussian field at early 
times, as required from the inflationary paradigm of mod¬ 
ern cosmology. Typical hierarchical amplitudes such as the 
time-honored one-point cumulants Sp (Bernardeau 1994) do 
match well on these scales those expected in a ACDM Uni¬ 
verse, showing consistency with PT. Early works have shown 
that the lognormal one-point probability distribution func¬ 
tion (PDF) is thus a good match to that observed in simu¬ 
lations (Bernardeau & Kofman 1995; Kayo et al. 2001, e.g.). 
More recently, the lognormal field was proved successful at 
modeling much more elaborated properties, such as the cou- 
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pling of very large wavelength perturbations to small scales 
Fourier modes (Takahashi et al. 2014; Carron & Szapudi 
2015). The ease with which it can be simulated makes it also 
a most practical and useful tool to simulate covariance ma¬ 
trices for cosmological data, and it is successful as a prior in 
bayesian matter field reconstruction from large scale struc¬ 
ture data. (Hilbert et al. 2011; Kitaura et al. 2010; Carron 
et al. 2014; Ata et al. 2015). 

One of the most intriguing property of the lognormal 
field S is that the full set of the in cosmology most tradi¬ 
tional statistics the A^-point functions (Peebles 1980; Fry 
1984; Bernardeau et al. 2002; Szapudi 2009) 

^jv(a;i, ■ • ■ ,xn) = {S{xi) ■ ■ ■ 5 {xn))^ (1) 

or their Fourier transform the polyspectra does not pro¬ 
vide a complete statistical description of the field, as already 
pointed out in (Coles & Jones 1991). This is generic behav¬ 
ior for skewed nonlinear fields in which the series expansion 
of the moment generating function is only asymptotic (Car¬ 
ron 2011), as is the case for the cosmological matter density 
field. We refer the reader to (Valageas 2002) for a very de¬ 
tailed study of the density field moment generating function; 
the main idea of how this can be is fairly simple: given the 
observed non-linear density J one may write schematically 
the moment generating function going back to the 

initial, Gaussian, linear density as 

(e-“) (2) 

If nonlinear growth J[(5n] more than compensates the Gaus¬ 
sian decay —(5| in the exponent, then the above expression 
diverges for any t < 0. This prevents the convergence of the 
series ~ /^' ^ however small, 

and thus the recovery of p{S) from the moments. 

For such fields, polynomials in the field do not form a 
complete basis set of functions and thus only part of the the¬ 
ory is probed by the hierarchy. It follows that there exists 
a variety of probability functionals that predict the same 
hierarchy of correlators, some of which given by (Carron & 
Neyrinck 2012), and information escapes the hierarchy. A 
deeper physical interpretation of this effect remains myste¬ 
rious however, and it is not entirely clear what it is that the 
hierarchy does and does not characterize. It is the purpose 
of this Letter to elucidate these questions in the case of the 
lognormal field. 

We show here that the hierarchy (in fact, a some¬ 
what larger set, defined below, allowing in effect the use 
of Laurent-series in place of only ordinary power series in 
the field) does measure the probability functional not of the 
Gaussian log-density field, but that of a different random 
field, the ’wrapped’ Gaussian field. The time evolution of 
both functionals is given by the same equation, the diffu¬ 
sion equation, with identical initial conditions. The former 
diffuses in Euclidean space keeping its Gaussian shape, the 
latter however on the compact torus, reaching homogeneity 
in the long-term, the information contained in the hierarchy 
diffusing to zero quite literally. Eventually, this will allow us, 
as our main result in Eq. 24, to characterize the hierarchy of 
moments as capturing only a specific part of the log-density 
field. In order to do this, we will only need derive the hier¬ 
archy information matrix and comparing it to that of a field 
different to the log-density. 


Gaussian probability functional and diffusion equation 
One way to characterize the Gaussian regime, thus, linear 
growth in the Universe, is as already advertised through 
a diffusion process. Let pc denote the Gaussian PDF for 
the fluctuation field 5 at any number d of different spatial 
positions. We use for simplicity as ’time’ coordinate t the 
cosmological growth factor . Using the vector notation 
5 = ((5(a:i), • ■ • ,5{xd)), pa can be explicitly written as 


PG{5,t) 


(27r)‘^/2.ydett^ 


(3) 


with covariance matrix = ^{xi — Xj) the linear two-point 
function of the field (at t = 1). The time evolution of pc is 
given by the diffusion equation 


dpG _ 1 ^ d'^PG 

dt ^ dSidSi ’ 


(4) 


with diffusion constants matrix The lognormal field pre¬ 
scription amounts to the replacement of the Gaussian den¬ 
sity fluctuation field J by a Gaussian log-density field A = 
ln(l -I- J). The most trivial benefit is that in contrast to the 
Gaussian case, the constraint J = p/p — 1 ^ —1 is ensured 
in a sensible way. Mass conservation (the condition (J) = 0) 
is also ensured by forcing a nonzero mean (A) = —(t^/ 2, 
where is the variance of the log-density field. We set 
however the mean of the A field to vanish in what follows, 
for notational convenience. This is irrelevant to our results, 
the interested reader can introduce the correct value in each 
step at the cost of extra ink. 


2 N-POINT CORRELATORS IN THE 
LOGNORMAL FIELD 

We work with an arbitrary number d of pixels with coor¬ 
dinates Xi. The two-point function of the Gaussian A field 
on these pixels is the finite matrix ^A,ij = f,A(xi — Xj). We 
do not write explicitly a time dependence, keeping in mind 
^A oc t. For lognormal field statistics, the following con¬ 
venient, compact equation holds true for a vector n with 
entries Ui taking arbitrary values (Carron & Neyrinck 2012) 
((^) = 0), 

= exp • n) . (5) 

The standard hierarchy of A-point moments is described by 
the set of such moments for which the entries Ui take all pos¬ 
sible non-negative integer values. We extend the hierarchy, 
allowing negative values as well: the vector n (’multiindex’ 
henceforth) runs in what follows over lattice points .The 
reason will be clear in a moment. The set of observables we 
consider is thus 

(1 + <5(a;i))’^A , with n G (6) 

i=i / 

Let further be for each pair of multiindices n, m the ma¬ 
trix element Mnm be defined through Mnm = mn+m. The 
matrix element Mnm — mnmm is the covariance between 
correlators rUn and mm. As a meaningful information mea¬ 
sure we consider Fisher information (Tegmark et al. 1997, 
e.g.), an information measure associated to arbitrary model 



n(l + J(a:0r) 
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parameters Q, /5, that quantifies how well these parameters 
may be inferred from the data under consideration. Our aim 
is to evaluate the information contained in the hierarchy (6) 
of observables, defined by 


T^hierarchy 

-^a/3 


r^—l-j d^rim 
da ^ Jnm Qp 

n,m in 


( 7 ) 


which is at most that contained in the entire (density or 
log-density) field given by 

,:,tota,l f jd ^daPG[A)dppG[A) 

The former reduces to the latter when the probability func¬ 
tional is uniquely defined by the moment hierarchy (see 
(Carron 2011) for a discussion). is in our case the in¬ 

formation of a Gaussian random field, given by a well-known 
formula (Vogeley & Szalay 1996; Tegmark et al. 1997) col¬ 
lecting the independent Fourier modes of the field; we will 
not need here its specific form however. 


2.1 Inverse cosmic covariance matrix 


We now obtain the inverse matrix M~^. This can be per¬ 
formed elegantly with this multindex notation as follows. 
We first note the simple transformation 


i (n-f m)^CA • (n-fm) = 

n^^A • n -f m^^A • m - i (n - m)^ ^a • (n - m) . 


We then observe that 
exp 


--(n-m) ^A-(n-m) 


rn-n—m 

Using these two relations in Eq. (5), we can write 


( 9 ) 


( 10 ) 


( 11 ) 


We then define the matrix Mnm = l/nin-m. This matrix is 
invariant under lattice translations n, m — >■ n + r, m-f r, and 
this allows its inversion below through Fourier transforma¬ 
tion. Writing relation (11) in terms of the matrices M and 
M, and inverting on both sides, results in 


[M-] 


1 1 
m2 m^ 



( 12 ) 


It remains to invert M. To that aim, we introduce its Fourier 
transform. Let 6 = {6i, ■■ ■ , 6d), and 


pw{6) = 




rrin 


(13) 


From (10), this series is absolutely well behaved and we need 
not worry about convergence issues. Inverting this relation 
gives the Fourier representation of M, 

Mnm = ^^= [ d‘^epw{0)e*^'^^-’^h (14) 

The region of integration is for all components 9i in [—7r,7r] 
(or any interval of length 27r), a short hand notation for 
which is the d-torus T'*. This provides us with the inverse 
matrix : 



d‘^6 


je-(n-in) 


pw{d) 


(15) 


A last intermediate result is required before the evaluation of 
the information matrix. Consider the impact of a parameter 
a on pw{9)- From its Fourier representation follows directly 


dpwjd) 

dot 



OlTlti 1 ptti-6 

da rrin 


(16) 


We now reach the main result of this Letter. From Eqs. 
(7),(12),(15) and (16) we obtain a remarkably simple ex¬ 
pression 


^hierarchy 

^a/3 



d'^e 


daPw{9)dfiPw{9) 

Pw{6) 


T-itotal pwr 


(17) 


It has exactly the form (8) of the total information content 
not of the Gaussian field A but of a new field 9, on the torus, 
obeying the PDF pw- 


2.2 The ’Wrapped’ Gaussian field 


We turn to the justification of our our interpretation of pw 
as a new random field PDF and its physical meaning. We 
derive another useful representation of pw, that makes its 
connexion to pg clearer. Written in full, Eq. (13) gives the 
following 

^ expj^-in^eA-n-ie-n^ (18) 

n in 


We transform the first exponential factor using the standard 
Gaussian integral formula 


a-WiA 




d'^A 




■A-\-iA-n 


gd (27r)‘'/2.^det^A 


(19) 


The sum over the multindices in Eq. (18) reduces now to 
the d-dimensional Dirac-comb 

1 ^ ^ 5^i2.n-{x-9)). (20) 

n in n in 


The Dirac delta allows us to undo the x-integrations, with 
the result 


Pw{9) 


1 

(27r)‘*/2.ydet ^a 



g-i(S-2,rn)'^«-l.(S-2,rn) 


( 21 ) 

This makes the positivity of pw obvious. On the other hand, 
representation (18) tells us that pw integrated on the torus 
gives unity. This justifies fully our interpretation of pw as 
a random field PDE. We note from this representation (21) 
that Pw reduces to the Gaussian field at early times, given 
by the term n = 0. 

Now, what is the physical interpretation of this new field 
? In fact, Pw satisfies exactly the same diffusion equation 
than Pg- Remember that ^a evolves proportional to t. From 
its representation (18) follows indeed 


dpw _ 1 ^ d^pw 


( 22 ) 


where ^a in this equation is again the two-point function at 
t = 1. From representation (21), initial conditions are the 
same, the only difference is topological: the field 9 diffuses 
on the torus. It is physically obvious that the topological 
constraint must change the longterm behavior of the field. 
In particular we expect pw to reach homogeneity, as soon 
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as, say, a a ~ tt. From representation (18) holds indeed that 
for any shape of two-point function holds at late times 

^ =^nd ^ 0. (23) 

This is a rigorous demonstration that the joint information 
of the hierarchy decays to zero. The critical time cta ~ rr is 
larger than a previous estimate a a ~ 0.4 one of us obtained 
considering one-point moments only (Carron 2011). This is 
due to the presence of the inverse powers in the hierarchy 
considered here. 


3 FINAL WORDS 

We are now in a position to answer the title’s question in 
the case of the lognormal field. A rigorous statement is as 
follows. In the log-density field decomposed at each point 
into a phase angle 0 in —tt and tt and an amplitude k (an 
integer) 

ln(l + (5(x)) = ^(a;)-I-fc(a;) 27r, (24) 

the hierarchy (6), including powers and inverse powers, mea¬ 
sures only the phase held 6. This may as well be formulated 
in terms of a lossy nonlinear transformation, with which we 
conclude. We consider two hypothetical analysts, facing data 
in form of the lognormal held on d pixels. The hrst extracts 
the hierarchy of correlators (6) and discards he original data. 
The second reconstructs the phase held the log-density in 
each pixel, 

ln(l + (5(xi)) —>• 6l(xi). (25) 

This can be pictured as wrapping the real line on the unit cir¬ 
cle. Following the rules of probability theory, this new data 
is described precisely by the wrapped Gaussian held pwifi), 
Eq. (21). This transformation makes no difference in the ini¬ 
tial conditions. Then, under-densities start to permeate in 
the most over-dense regions, and vice-versa. Eventually, any 
memory of the initial held is lost, and the data can perfectly 
be described by white noise. Yet, both analysts have the 
same information matrix, and are going to reach identical 
conclusions when performing inference on parameters from 
their data. 

Sourcing this effect is the rapid growth of the moments 
across the hierarchy, as the density held evolves into the 
cosmic web. Clearly, the matter density is not exactly log¬ 
normal, but, for the reasons exposed in the introduction, it 
must also be affected by this phenomenon. Trying to under¬ 
stand it for a more realistic model might well be a rather 
formidable task. This is left for future investigations. 
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